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The optical conductivity of graphene strained uniaxially is studied within the Kubo-Greenwood 
formaUsm. Focusing on inter-band absorption, we analyze and quantify the breakdown of universal 
transparency in the visible region of the spectrum, and analytically characterize the transparency 
as a function of strain and polarization. Measuring transmittance as a function of incident po- 
larization directly reflects the magnitude and direction of strain. Moreover, direction-dependent 
selection rules permit identification of the lattice orientation by monitoring the van-Hove transi- 
tions. These photoelastic effects in graphene can be explored towards atomically thin, broadband 
optical elements. 
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I. INTRODUCTION 

Transparent flexible electronics is currently a much 
sought technology, with applications that can range from 
foldable displays and electronic paper, to transparent so- 
lar cells. Graphene, on account of its ultimate thick- 
ness, large transparency [1], high mechanical resilience 
under strong stress/bending cycling, and excellent elec- 
tronic mobility [2] , has been promptly ranked among the 
best placed materials to achieve those technologies [3, 4]. 
Such goals require a thorough understanding of the inter- 
play between the different aspects that will unavoidably 
be present in such devices, namely how sensitive the di- 
electric and optical properties of graphene are to gating 
and straining. At the same time, broadband optical el- 
ements that can be scaled down to the nanoscale, and 
easily integrated into photonic/photoelectronic circuits, 
are equally appealing scenarios in current nanotechnol- 

ogy- 

The optical absorption response of graphene has been 
recently given thorough attention on both the experi- 
mental [1, 5-8] and theoretical fronts [9-17]. One of 
the distinguishing features of undoped (or lightly doped) 
graphene arises from the constancy of its transparency, 
T(w), which is controlled by electron- hole excitation pro- 
cesses, and universal: T{uj) w 1 — 7ra (a ~ 1/137 being 
the fine structure constant) [1]. This universality is a 
consequence of the particle hole symmetry of graphene's 
spectrum, combined with the cancellation of the fre- 
quency dependencies of the matrix element and vanish- 
ing density of states. It holds throughout a broad spec- 
tral range comprising the frequencies between the Fermi 
energy, /z, and the vicinity of the van Hove singularity 
(VHS) at the M point in the Brillouin zone (BZ). In un- 
doped graphene this covers a band spanning the visible 
region, down to the far IR. 

Here we analyze and quantify how the optical ab- 
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FIG. 1: Lattice orientation in direct and reciprocal space, 
overlaid with several quantities used in the text. The bottom 
illustration shows the relative orientation of the lattice with 
respect to a general tension direction 9. 



sorption associated with direct optical transitions is af- 
fected by strain. The strain-induced anisotropy leaves 
a clear signature in the optical response of the system, 
modulating its transmittance, reflectance and absorp- 
tion, while simultaneously rotating the polarization of 
incoming light. Tailoring of graphene systems on the 
basis of such properties extends the concept of strain 
engineering in graphene from electronic structure and 
transport [18, 19] to the optical domain as well. Such 
effects are described next. Related work on effects of 
strain in graphene's optical response has been reported 
in Refs. [20-22]. 
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II. STRAIN-INDUCED ANISOTROPY 

We concentrate on uniform uniaxial deformations of 
free graphene, which provide the highest degree of 
anisotropy for given amount of strain, e. Within a tight- 
binding approximation for the 7r-band subsystem, strain 
impacts electronic motion via the modification of the 
Slatcr-Koster parameters = t(Si) in the Hamiltonian 



= h'^VpXq.T]'^ .q, where hvp = 3ta/2. Diagonalization 
of T]^ yields the principal velocities v^. = Vprj^., with 



+ ^2 ~^ ^3 
3t2 



(1) and the principal directions: 



R,S 



Here R denotes a position in the Bravais lattice; ^1,2.3 
connects site R to its neighbors; a{R) and b{R) are 
the field operators in sublattices A and B. In what fol- 
lows, we characterize the interplay between strain and 
electronic structure in the same framework described 
in detail in ref. 23. In summary, this entails the as- 
sumption that hopping amplitudes vary with distance as 
t{r) ~ ie-^('-"\ with t w 2.7 eV, a = 1.42 A, Aa w 3 - 4 
[21, 23], and we disregard bond bending effects, since 
they are not significant in this effective description [24]. 
Bond deformations are uniform and given to linear or- 
der in terms of the strain tensor Sij: 6 — {1 + e) ■ 5'^. 
This approximation was shown to describe with good 
accuracy both the threshold deformation for the Lif- 
shitz and metal-insulator transition at large deformations 
[23, 25, 26], and the behavior of ti{5) when compared to 
ab-initio calculations [24]. Since the hexagonal lattice is 
elastically isotropic, Sij can be fully parametrized by the 
amount of uniaxial strain, e, its direction with respect 
to a zig-zag (ZZ) direction, and Poisson's ratio, i' « 0.16 
[23]. Our Cartesian directions are such that a ZZ edge 
coincides with the x axis (fig. 1). Tension along 6 = 
[ZZ] and 6 = 7t/2 [armchair (.A^)] defines representative 
directions that will be recalled frequently. We shall dis- 
regard the strain-induced modification of the reciprocal 
lattice, as it is not relevant for the optical absorption. 
Since here we are not interested in large deformations, 
we further approximate U = t{Si) « i — tX{Si — a). 

The low energy Hamiltonian appropriate for optical 
processes below the UV can be derived conventionally 
[16], by expanding (1) around the shifted Dirac points: 
k = kjj + q {q <ti ko)- For arbitrary ti the nonequivalent 
Dirac points lie at ko.^ = ({A ~ B, A + B), with A = 



V3a ''^^^ 2*2*3 ' ^ 3a 2ti*2 

is the valley index, labeling the two nonequivalent Dirac 
points. In their vicinity, the electronic dispersion is given 
by 



-1 ti 



B 



and C = ±1 



^1 



3\/3 
3 



q^qya [t^ ~ t^) 



+ -qia'{2ti~ti + 2ti)+Oiq'). (2) 

The Fermi surface is consequently an ellipse, with prin- 
cipal axes that will be rotated with respect to Oxy in 
the general situation, eq. (2) can be cast compactly as 



tan(y9± 



+2 _ +2 

(.3 



Viivlt^-tlY 



(4) 



These directions define the slow/longitudinal (—,0 ^^id 
fast /transverse {+,t) directions of strained graphene, in- 
sofar as the direction-dependent Fermi velocity in the 
elliptical Fermi surface is smallest along the one, and 
largest along the other. In the principal coordinate sys- 
tem, the effective Hamiltonian of valley C — ^ reads 



(5) 



where ri_2 are Pauli matrices acting on the {A, B) sub- 
lattice space. Coupling to a light wave described by 
the physical vector potential A(t) = Aoexp(— icji) is 
achieved by the minimal substitution p — > p + eA in 
eq. (5). The frequency-dependent conductivity is ex- 
tracted from the linear response to Aq. 

For future reference, the anisotropy parameters will be 
expressed directly in terms of the longitudinal deforma- 
tion e, to first order. For general tension along 9 with 
respect to Ox we have 



77t ~ 1 -I- aXve, 77; ~ 1 — aXs, 



5k 



CA£^-4^(cos26',-sin26'), 



(6a) 
(6b) 



tan ifii 
tan ipt 



l-aA£^-^(l + 2cos46l) 
8 

l + aAe^-^(l + 2cos46l) 



tan^, (6c) 



tan (6* 



(6d) 



As intuition dictates, the slow/longitudinal axis is co- 
incident with the tension axis {(p w 6). From (5) the 
strain-induced correction to the density of states (DOS) 
in the vicinity of the Dirac points follows immediately: 
p{E) ~ p^^° (E) / (riiriy) , where the isotropic DOS reads 
p^''°{E) = 2\E\/{Trh'^vl,). Notice that, even though the 
effective longitudinal (transverse) velocity increases (de- 
creases), the net effect in the DOS is always a slope en- 
hancement, because rji'qt ~ 1 — aA(l — i')s < 1. This 
means that strain modifies the Fermi energy and/or elec- 
tron density. 
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III. CONDUCTIVITY OF STRAINED 
GRAPHENE 

Kubo's form of the frequency dependent conductivity 
reads 



^2 ykXX' ^,kX' X 



(7) 



where gs = 2 is the spin degeneracy, A the total area, 
A, A' = ±, and v^^^ = (fc, A|wct|fe, A') are the matrix 
elements of the velocity operator Vq. = ih^^[H,Xa] in 
the momentum eigenbasis. For the general tight-binding 
Hamiltonian (1), these matrix elements read explicitly 



A 



of the Dirac approximation, and for 2/i < ^ T, be- 
ing only renormalized by the anisotropy factors 7y± . The 
degree of anisotropy is controlled by rji/rjt, which is a 
sensible result because crii(w) cx; |Uj ' '~ and the ratio 
reflects the quotient between the Fermi velocities along 
the principal strain directions. From Eqs. (10) and (11), 
and recalling the expressions in eq. (6a), we can express 
the strain induced corrections to the full isotropic con- 
ductivity, cr(w) = o''('^) + *o'"('^)i in linear order in the 
deformation as 



\uj) X (1 T2|5fcD|a) 



(12) 



where a^^°{uj) represents the full frequency-dependent 
conductivity in the absence of strain. For example, ten- 
sion along ZZ {9 = 0) yields a decrease in Uxx, and an 



k,\\v\k ,A ) = — — S^^k'Sw ^Q^a sin [fc.ricj — 0^] increase of the same magnitude for ayy. Substitution 



iX 



Sk.k' (1 - ^ ta6a COS [k.Tla - 9k] , (8) 



with 9k = arg[^j exp(jfc.ni)] . Contributions with 
A = A' are associated with mtra-band transitions, and 
A 7^ A' with mter-band, direct transitions. In this form 
the conductivity can be directly calculated at the tight- 
binding level. In order to proceed fully analytically we 
work in the Dirac approximation (5), and consider a clean 
system at zero temperature, thus retaining only the inter- 
band contribution. In the principal coordinate system de- 
fined by the slow/longitudinal and fast /transverse axes 
the relevant matrix elements become: 



fc,A,-A|2 



fc,A,-A|2 



fc.A,-A fc.-A,A 



• vlrjf sin^ 9kr. , (9a) 
' vl 1]^ cos^ 9ko , (9b) 
;^42^sin(2(?fc„). (9c) 



The anisotropy is explicitly encoded both in the param- 
eters r?±, and in the phase Oki^ ■ 

From here and eq. (7) it is straightforward to obtain 
the frequency dependent conductivity CTq,^ = a'^^ + ia'^ 
Its real part reads 



/(-¥-A*)-/(¥ -A^) (10) 



for the longitudinal conductivity, and (j'tti.'^) — 
{rjt / riif' a'liiiu) for the transverse component. The 
isotropic (and universal) value is ctq — /{Ah), and f{x) 
represents the Fermi-Dirac occupation function. The 
imaginary part, a"{Lo), reads 



— X X lo^ 



21/^1 



■ w 



2\^i\ 



(11) 



when T = 0, and cf^^{uj) — {"nt / rjiY a'^. The off-diagonal 
components au are zero, as one expects from symmetry 
and the absence of magnetic fields. This result shows 
that cr^j(w) is still constant within the region of validity 



vv 

of the material parameters applicable to free-standing 
graphene in (12) and (6b) yields an anisotropy factor 
{(Jt — (Ji) / (T^^° ~ A\5kD\a ^ Se. This is sufficiently marked 
to be detectable by conventional optical means, like ab- 
sorption in the visible or IR, which we discuss below. 



IV. AB-INITIO OPTICAL CONDUCTIVITY 

Given the approximations used, it is legitimate to ques- 
tion the validity of the analytical corrections written in 
eq. (12). To that end, we have extracted the optical 
conductivity of uniaxially strained graphene from first- 
principles as well. Density Functional Theory (DFT) cal- 
culations were performed with the ab-initio spin-density 
functional code AIMPRO [27]. We used the GGA in the 
scheme of Perdew, Burke, and Ernzerhof [28] . Core states 
were accounted for by dual-space separable pseudopoten- 
tials [29]. The valence states are expanded over a set of s, 
p, and d-like localized, atom-centered Gaussian functions. 
The BZ was sampled according to the scheme proposed 
by Monkhorst-Pack [30]. We used a convergence-tested 
grid of 20 X 20 X x 1 points for the self-consistent calcu- 
lations. In equilibrium we obtained the optimized lattice 
parameter a = 1.42 A. 

Uniaxial strain was applied with relaxation as de- 
scribed earlier in Ref. 24. For each strained configura- 
tion we extracted the dielectric function within the long- 
wavelength dipole approximation [31], and from it the 
optical conductivity. The results for (t^^(w) and a'yy{uj) 
so obtained are shown in fig. 2, for representative uniaxial 
strains applied along ZZ {9 — 0) and AA {9 — 7r/2). 

The main panel of fig. 2(a) shows the renormalization 
of the real part of axx, which corresponds to an for ZZ 
and att for AA. It is clear that for visible and IR fre- 
quencies the conductivity remains roughly constant in 
uj, but its magnitude depends on the amount of strain. 
The strain dependence is shown in detail in the inset, 
at = 1 eV. The perfect linearity in e of the ab-initio 
results up to at least e ~ 10% shows that the analyt- 
ical result of eq. (12) is indeed quite dependable for a 
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FIG. 2: Real part of the optical conductivity a:cx (a) and Oyy (b), calculated ab-imtio for graphene under uniaxial strain 
along ZZ (6 = 0) and AA {0 — it/ 2). The sharp absorption peaks are associated with transitions resonant with the van Hove 
singularities. In panel (a) the peak is resonant with the transition at Mi, 3 in the BZ, whereas in panel (b) the double peak 
structure derives from the splitting between transitions at Mi_3 and M2. Both insets show the variation of a'^^ with strain at 
Lj — 1 eV, as obtained ab-mitio. The straight lines are linear fits, which yield slopes of —2.45 {ZZ) and 2.64 (.4.4) for the inset 
of panel (a), and 3.22 (ZZ) and —2.08 (.4.4) for the inset of panel (b). 



wide range of stretching. We must point out, however, 
that the linear slopes quoted in fig. 2 cannot immedi- 
ately be used to extract the bandstructure parameter A 
using eq. (12). That is because the ab-initio calculation 
naturally includes the relaxation and deformation of the 
lattice, which, as advanced in the very beginning, we 
chose to ignore, not to encumber the discussion by more 
complicated expressions which do not change the main 
results. 

Taking the lattice deformation into account leads to a 
correction of the expansion (2) of the dispersion around 
the Dirac points. This leads to a renormalization of en- 
tries in the matrix 77^ . It can be shown straightforwardly 
that such changes amount to replacing Aa — > Aa — 1 in 
t]^ , and, consequently, in every ensuing result. In this 
way, the ab-initio slopes quoted in fig. 2 correspond to 
A ~ 2.42 A"^ (ZZ) and A ~ 2.56 A"^ (AA) for the inset 
of panel (a), and A ~ 2.97 A"^ (ZZ) and A ~ 2.17A"^ 
{AA) for the inset of panel (b) . This rather nicely agrees 
with the estimates quoted earlier that put Aa ^ 3 — 4. 



V. VAN HOVE SINGULARITIES 

In the UV band, cr{uj) is dominated by direct tran- 
sitions between the AI points in the BZ, which coin- 
cide with VHS in the electronic dispersion. From the 
point of view of (1), these transitions occur at momenta 
k coinciding precisely with Mi^a — (±7r/-\/3, tt/S) and 
M2 = (27r/3, 0) (see fig. 1 for our convention regard- 
ing the M points). The corresponding resonant frequen- 
cies are therefore given by ujMi — 2\E{k — Mi)\, and 
read ujmi = 2|ti -\- 12 - iai, WAf2 = 2|ii - t2 + hi, and 



— 2\ ^ ti + t2 + t^l. In the most general situation 
ti t2 ts the VHS split. Linearly expanding the hop- 
pings in the magnitude of uniaxial strain, and defining 
ujMi ~ 2i -I- taXeAuJMi, such splitting acquires the ex- 
plicit form 

AwA/1,3 ~ I' - 1 + (1 + i^)(cos26l± \/3sin26'), (13a) 
AcjA/2 ~ I' - 1 - 2(1 + t/)cos26'. (13b) 

But, as is obvious from the ab-initio results in fig. 2, this 
splitting is not always visible in the absorption spectrum. 
That is because the velocity matrix elements (8) impose 
a modulation of the strength associated with these tran- 
sitions. In order to extract this effect we need to abandon 
the Dirac approximation (5,9), and work with the full- 
tight binding bandstructure and matrix elements. In (8) 
take, for example, the matrix elements at A; = M2 , which 
read 

K''^l' = |J(^i-*3f, (14a) 
K'\'-'^iti + 2t2+ts)^ (14b) 
= + 2t2 + h){t, - t,). (14c) 

When expanded in powers of strain, will be zero 

to linear order, irrespective of the strain direction. This 
suppresses the M2 singularity in cr'(i^) leaving only the 
peaks related to Mi^s, whose associated matrix elements 
are finite. Since the geometry chosen in fig. 2 consid- 
ers only 9 = and 7r/2, Mi and M3 are still degener- 
ate, as per (13), and thus only one peak should survive 
in (T^j,, precisely as seen ab-initio in fig. 2(a). On the 
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that would be 



tan^T 



tan ( 



Kil-4\SkD\a— , (15) 



FIG. 3: Illustration of the behavior of the transmittance as 
a function of incident polarization measured in the labora- 
tory frame. According to (16), the phase shift 6 defines the 
direction of strain, and the amplitude its magnitude. 



other hand, there is no such selection rule arising from 
the velocity matrix elements when computing a'yy, which 
allows the splitting between M2 and Mi_3 to be observed, 
as fig. 2(b) clearly demonstrates. Eqs. (13) also explain 
why in fig. 2(a) the strain- induced shift of the peak is 
less pronounced for tension along ZZ (9 — 0), than AA 
{9 = n/2): Aujmi{ZZ) = v Awmi(^-A). And similarly, 
Eqs. (14) account for the different relative intensity of 
the peaks associated with Mi and Mi^a in fig. 2(b). 

It is important to stress here that these results are not 
tied to a particular coordinate system. We can easily ex- 
press Eqs. (14) in any coordinate system, rotated by Lp 
with respect to Ox'\i^ and conclude that the longitudinal 
matrix element | v^, ' | ^ always vanishes when coincides 
with a ZZ direction. For example, vanishes for 

</? = 0, l^i^^P vanishes for = 7r/3, and 1^*/^^ vanishes 
for Lp — l-K I'i. Suppression of a van Hove peak in the lon- 
gitudinal conductivity therefore singles out one of the ZZ 
directions of the lattice. This is consistent and explains 
the numerical calculations of Ref. 21, and has an im- 
portant consequence: by monitoring the splitting in the 
absorption peaks associated with VHS, and by inspecting 
the selection rule just described, one can measure simul- 
taneously: the magnitude of strain, its direction, and the 
direction of the underlying lattice with respect to the labo- 
ratory coordinate system. This provides an optical means 
to perform the same measurements that have been made 
by exploring the splitting of the G peak in the Raman 
spectrum of strained graphene [32, 33]. 



VI. TRANSPARENCY AND DICHROISM 

The asymmetry (12) in the conductivity tensor will re- 
sult in a certain degree of dichroisni, as the absorbance 
of linearly polarized light will depend on the polarization 
direction with respect to the slow/fast axes. Treating 
graphene as a 2D conducting sheet, and solving the as- 
sociated Fresnel equations, we can extract the degree of 
polarization rotation for normal incidence on graphene 
in vacuum. In the visible and IR where eq. (12) stands 



where (j)T.i are the transmitted and incident polarizations 
measured with respect to the slow/longitudinal axis, and 
cfj,oao = 7ra ~ 0.02. Likewise, the transmittance for 
linear polarization becomes 



T ~ 1 — 7ra [1 — 2 1 I a cos 2(/)/] 



(16) 



An important consequence of this periodic modulation is 
that it allows a direct determination of both the strain di- 
rection and its magnitude, as follows. In some laboratory 
coordinates (16) is transformed by making — ipi. 

The amplitude of T as a function of polarization direc- 
tion determines the amount of strain, while the phase 
shift (fi fixes the direction. This is illustrated in fig. 3. 
The corrections to both polarization and transmittance 
are weighted by ira, and will be necessarily small. But, 
one the one hand, the modulation amplitude is roughly 
^ Seira (0.16£) and transmittance can be routinely mea- 
sured within 0.1% of precision. On the other hand, the 
transmittance of multilayer graphene is, to a great ex- 
tent, cumulative [1, 7]. This implies that the same re- 
sults apply for multilayer graphene, provided one replaces 
ira — ^ Nna in (16), where N accounts for the number 
of graphene layers. The effect is thus naturally enhanced 
in multilayers, as it is in the vicinity of the VHS (or any 
resonance, for that matter). Another interesting appli- 
cation is that, if strain can be controlled with precision, 
an expression like (16) allows, by means of a simple op- 
tical experiment, the measurement of the bandstructure 
parameter A — d\ogt/dr, whose knowledge is crucial for 
the characterization of all strain-induced effects on the 
bandstructure. 



VII. DISCUSSION 

The photoelastic effect in undoped graphene has been 
quantified, and shown to possess features that might be 
appealing in the development of atomically thin optical 
elements. One of such characteristics is the frequency 
independent response in a very large frequency range, 
which remains valid when the system is anisotropically 
strained. This constancy and predictability is a rele- 
vant feature for broadband applications. The degree of 
anisotropy induced by strain is determined by how much 
the Dirac point is displaced from its position at K/K' 
in the BZ: \5kD\ (6b). This opens several possibilities, 
such as: monitoring optical absorption as a function of 
strain to characterize the band parameters of graphene; 
or monitoring the transmittance as a function of incident 
polarization in order to measure the magnitude and di- 
rection of strain in graphene devices. This last example 
could find applications in completely passive, transpar- 
ent, strain sensors. 
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Even if the magnitude of the photoelastic effect in 
graphene at visible or IR frequencies is Umited by its 
small natural absorption of 7ra ^ 2%, it is nonetheless 
significant for an atomically thin membrane. Additional 
versatility is provided by the fact that the effect can be 
naturally amplified by stacking multilayers, or that the 
optical absorption can be radically affected by electronic 
doping as well, on account of Pauli blocking [5, 9]. 

A determination of the lattice orientation cannot be 
made from the optical absorption at low energies, sim- 
ply due to the isotropy of the Dirac dispersion. By con- 
trast, the electronic dispersion at energies close to the 
VHS fully reflects the symmetry of the underlying lat- 
tice. Therefore, by measuring the optical response at 
frequencies resonant with the van Hove transitions, and 
analyzing the strain-induced splitting and selection rules 
of the absorption peaks, one can extract the amount of 
strain, its direction, and the lattice orientation. 

Another consequence of the strain-induced tunability 
of the VHS is that it can have a significant import in cur- 
rent efforts to elevate the Fermi level of graphene up to 
the VHS [34-36]. The ability to achieve this, and thereby 
dramatically increase the electronic DOS, is expected 
to facilitate many-body instabilities and the establish- 
ment of correlated phases, such as superconductivity, or 
charge/spin density waves [37, 38]. Strain engineering of 



graphene can, in this respect, work as a facilitator and 
provide added tunability. For example, from (13) it fol- 
lows that the energy of the VHS can be reduced by uniax- 
ial strain to values as low as E/t ~ 1— (3-|-i^)aA£/2. Using 
Aa ^ 3 for estimate purposes we can write E /t c^i 1 — 5e, 
so that the reduction for 10% strain can be as much as 
50% in the position of the VHS. On top of this we would 
need to include excitonic corrections that are known to 
renormalize the VHS further down in energy [39], and 
which we neglect in our treatment. 

Finally, the photoelastic effect is the basis of many op- 
tical and mechano-optical devices and applications at the 
macro-scale. The characteristics of graphene in this re- 
spect, which we just surfaced, might provide a valuable 
route towards the downscaling of those concepts to the 
realm of atomically thin optical elements, and their ap- 
plication at the nanoscale. 
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